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Design of Lecture 

It is intended to present incompressible flow CFD from the real-world 
applications point of view. It will be organized in three parts: 

Part 1 : Basics 

Formulations, solution methods, and some historical notes 
Part 2 : Applications 

If the value of CFD tools is viewed from engineering practices, there are 
issues beyond flow solution methods. The process is discussed using 
following two examples: 

- Turbopump 

- Biofluid 


Part 3 : Discussion 

Discussion/seminar on "Current Topics in CFD" 

Much discussion has been on-going whether CFD is a "mature technology 
A personal view on the state-of-the-art in CFD will be discussed 




Outline of Presentation 


• Introduction 

• Review of Solution Approach 

• Artificial Compressibility Method 

- Steady-state formulation 

- Time- accurate formulation 

• Pressure Projection Method 

• Solver/Tool development 

• Problem Solving Procedures 

- Historical example 

- Some building block studies 

- Parallel Implementation 

- Turbopump flow 

• Biofluid 

• Discussion/seminar on "Current Topics in CFD" 


Introduction 


• Background 

- Flow devices tend to be compact and highly ef f icient 

- Experimental approach can be expensive, 
and model tests need to be extrapolated 

- CFD may offer an alternative to reduce cost/time for development 


Scope of Presentation 

- CFD is viewed as an engineering tool 

- Will review solution methods for incompressible N-S equations with 
and emphasis on 3-D applications 

- Will discuss numerical/physical characteristics of primitive variable 
approach 

- Solution processes are discussed using real examples 




Introduction 


• Role of CFD 
Post 

Functionality was the most important aspect: 
Primary concern was "performace* prediction 

Current interest 


For engineering applications: 

"Cost" and "safety" are big issues -CFD is being used to reduce design 
cycle time and to expand safe operating envelope 

- CFD is being used routinely in aerodynamic design of commercial aircraft 
=> Return on CFD research is viewed as only incremental 

- Design of aircraft engine components is validated by CFD 

=> CFD- based design of an entire engine is yet to be realized 

- Virtual f light/maneuver and systems analysis require generation of large data sets 
=> Need CFD+IT : high-fidelity CFD solutions are still expensive 

For research (e.g. flow physics study or new enabling technology): 
"Idea" is the key, and can be explored via ‘big* problem 


Viscous Incompressible Flow 


# Formulation 

Can be viewed as a limiting case of compressible flow where the flow 
speed is insignif icant compared to the speed of sound 

- Low speed aerodynamics => Preconditioned compressible N-5 eq. 

Or truly incompressible 

- Hydrodynamics 


• Some Examples 

- Components of liquid rocket engine 

- Hydrodynamics (Submarines, propellers, ...) 

- Ground vehicles (automobile aerodynamics, internal flows...) 

- Biofluid problems (hemodynamics, artif icial heart, lung, ...) 

- Some Earth Science problems (with variable density) 





Solution Methods for Incompressible N-S Equations 


• Time- integration scheme: 

Primitive Variable Methods 

Based on Compressible Flow Algorithm 

- Artificial Compressibility Method (Chorin. 1967; Temam, 1977) 

ADI Scheme. FD (Centrakdiss) (BeamA Warming. 1978; Briley-McDonald, 1977) 

LU-5GS. FV (Centrakdiss) (Yoon 1987...) 

Line Relaxation. FD (Upwind) (... MacCormack, 1985) 

GMRES, FD (Upwind) 

Based on Pressure Projection 

- MAC (Harlow and Welch. 1965) 

- Fractional Step Method (Chorin. 1968; Yanencko, 1971; Marchuck, 1975...-) 

- SIMPLE type Pressure Iteration (Caretto et al„ 1972; Patanka A Spalding. 1972-.) 
Derived Variable Methods 

Vorticity- Velocity (Fasel. 1976; Dennis et al„ 1979; Hafez et a!., 1988 
Stream function-vorticity 

• Spatial discretization scheme: 

Finite Difference 

Finite Volume 


Governing Equations 


Incompressible Navier-Stokes Equations 
du. 

dt ox. ox. oxj 
Where s = source term 

du. du : 

r.. = 2 vS.-R..=vSl+^J- 

IJ IJ IJ ox . OX. IJ 

l J l) 


Eddy viscosity model for Rg 


R.=\r,.S..-2v.S.. z..=jfv+v t 'ti..=2v 7 S. 

ij 3 kk ij * ij ij \ tJij T ij 




GOVERNING EQUATIONS 
O Coordinate Transformation 

£ = £(*> y, *,*)>*? = r)(x, y, z, t), C = C(*> y» z > *) 
• Jacobian of the transformation 

H ««.->,<) _ i’ i> 

J — ■ det-r-; r — T)x 1}y *?z 

S(T,y,z) ^ (t (i 


• Metric terms are 


Vv z C - y< z v 

X^Zjf *^ 17 ^ 

x vy< - x <yv 


etc... 


GOVERNING EQUATIONS 


© Governing Equations 

^ * & li. i \ 


= - to - fa - to - 


where 


u = — v 

w 


1 r^p+«t/ 

e = 7 f»P + t»U 
47 £zP + wU 

U = Zt+txU + tyV+Zz w etc ‘ 


GOVERNING EQUATIONS 


© Viscous Terms 



d 

dxj 


uS, 





u x + U x 
V X + tiy 
w x 4 - u t 



Uy +V X 
Vy + Vy 
Wy + V z 


d 



u z + w x 

V Z +Wy 

w z + W z 


For constant v, the second terms in each bracket sums up 
to be zero 


GOVERNING EQUATIONS 


© Viscous Terms (cont’d) 

• For variable viscosity 



• For constant viscosity in orthogonal coordinates 



Artificial (or Pseudo-) Compressibility Method 


• Physical characteristics 

• Governing Equations in generalized Coordinates 

• Steady-State Formulation 
ADI Scheme 
LU-SGS Scheme 
Numerical Dissipation 

Artrf icial compressibility Parameter 
Boundary Conditions 
Geometry Effects, etc. 

• Time- Accurate Formulation 


Artificial (or Pseudo-) Compressibility Method 


• Artif icial Compressibility Formulation 
Continuity equation is modified to 



- Introduces hyperbolic behavior into pressure field. 

Speed of pressure wave depends on the artificial compressibility 
parameter, p. 

- The equations are to be marched in a time like fashion until 
the divergence of velocity converges to zero. 

=» Relaxes incompressibility requirement. 

- Time variable during this process does not represent physical 
time step. 





© Artificial Compressibility Relation 


l dp duj _ 

/? dr d£i 

• The equations are to be marched in a time like fashion until 
the divergence of velocity converges to zero. 

• The time variable for this process no longer represents phys- 
ical time. 

• In the momentum equations t is replaced with r , which can 
be thought of as a pseudo-time or iteration parameter. 




rU(i/-6)i 

£ 6 * p + uU 

J iyP + vU 

-tzP + V)U. 


where 
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NUMERICAL ALGORITHM 

© Time advancing 
• Linearization 

E n+ 1 = E n + A n (D n+1 - D n ) + 0 (Ar 2 ) 

where the Jacobian matrices are 

0 Li/3 L 2 P L$fi 

L\ Q + L\u L211 L3U 

I12 L\v Q + Ljv L$v 

-L3 h\ W I/2W Q "I" L3IO . 

Q = Zro + Li« + L2V + L^w 

Lo = te)«. £1 = £2 = (fi)„ £3 = ((,), 

(i = (f.1> or<)for(i,B, or 6) 


* 1 
i = w = i 


7T 


NUMERICAL ALGORITHM 


© Time advancing 
• Trapezoidal rule 


Z > n+1 =D n + 



+ 0 (At 3 ) 


Then the governing equation becomes 


D n+1 + =fj[6t(E - E v ) n+1 + Sr,(F - F v ) n+1 + ~ <X) n+1 ] 

= £>” - - £.)” + «,(# - F„) n + <<(6 - <?„)"] 
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NUMERICAL ALGORITHM 
O Governing equation in delta form 


{/ + \ J n+1 - Tx) + 6„{B n - r 2 ) + 6<(C n - r,)] | (£> n+1 - D n ) 
= -A rJ n+l [s ( (E - E,) n + 6 V (F - F,) n + 6 ( {6 - G.)"] 

+ (^ _1 ) Dn m - RHS 


O Approximate Factorization 


where 


L^L v L((D n+1 - D n ) = RHS 

b* 

11 

/ + T- JW+1 ^ {i "" 7l) ] 



L< = 

/+ ^ljn+^(C»_ 73 )] 


NUMERICAL ALGORITHM 
O Approximate Factorization 

• Second-order central differencing : 

- Block tridiagonal inversion 

(L„)AD = RHS 
(L<)AD = A D 

(L c )AD n+l = A D. 

! 

- Diagonal form is also available \ 

• Factorization error : cross-product terms 

h* [S^ASrjB + 6 n BS c C + S^C 6 (_A\ AG + 0 (h a ) 

where ’ „ ' ' 

A = A n - 7i, B = B n - 72, C — C n — 73, h=^J n+1 ■ * 
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NUMERICAL DISSIPATION (SMOOTHING) 
O Higher order smoothing terms are added for stability 


L i L v L((D n+1 - D n ) = RHS - €.[(V 4 A { ) a + (V,A„) 2 + (V ( A ( ) 2 ]Z?" 
where 


Li = 

/ + T JB+l ^ (in " 7i) + c ‘ v<A< ] 


1 + T jn+l6 v(B n ~ 7a) + <iV„A„] 

L< = 

/ + T jn+ltf<(d "" 73) + eiV<A< ] 


IMPLICIT vs EXPLICIT DISSIPATION 

O 1-D continuity equation : 5 ft + = 0 

O AF algorithm in delta form 

(assuming divergence free velocity at time level n) 

1 - ej V t A< (p" +l -p n j = -e e (VfA { )V 

O Discrete Fourier expansion of p 

P=£p(*)e ifc * 

n 

where 

k = — — n — wave number 

NA£ 

n = —N/2, ...0, 1, ..(iV/2 — 1) 

N = numiero/ mesh points 
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IMPLICIT vs EXPLICIT DISSIPATION (cont’d) 

Q This leads to [1 - «<*'] (p" +1 - p n ) = -e t (k') 3 p n 

where k' = -2 + 2cos(k) and ( k ') 3 = 6 - 8 cos(k) + 2cos(2k) 

O Amplification factor 

p n+1 [1 - t# - e t {k') 3 ] 
p” ~ [1 -**'] 

O To damp out numerical fluctuation 

kl < 1 

2e e < 


TT 


EXPLICIT DISSIPATION ON ACCURACY 

0 1-D continuity in AF form 

[1 - cjV, A,] (p n+l - p") = -f3hS x u - €«(V x A x )V 


O At steady state, this leads to 

6 Z u 


te(V,A x )V 

0h 


• Mesh refinement usually does not help 

• This error may be reduced (in a crude way) by 


where (e e ) 0 is the value from t = 0 to t = t Q 


Pressure Wave Propagation vs Viscous Effect 


© Locally linearized momentum equation 


^+2u-^ 

dt 2 atdx 

d 2 U d 2 U 

dt 2 dtdx 

These may be written as 


_ 3 &P 

^dx 2 

dx 2 


/a ^ r tp 

8t w 

~af 


' d_ 
dt 


+ (u + c) 


a_ 

dx 


d , 

m + {n 



( /38t w /8x\ 
y —dr w fdt J 


© Characteristic equation for linear waves 
(without shear stress terms) 


d , .8 

ai + (u+c % 


8 , 
* +(u 


C ) 


d_ 

dx 



= 0 
= 0 


» 

Pressure Wave Propagation vs Viscous Effect 


© One-dimensional characteristic equation (without viscous term) 


du 1 dp 
dt + (u ± c) dt 


du 1 dp 
dx + (ti ± c) dx 


(u ± c) = 0 


© Pseudo speed of sound 

c = \fu 2 + 


© Pseudo Mach number 



c 


u 

y/ti? +~P 


< 1 


=> M is always less than 1 for all /? > 0 





Criterion for (5 


© Two Time Scale 


• Time required for upstream propagating wave to travel a 
distance, x r 

XL 

tl = 


c — u 


• Time scale for viscous effect to spread a distance, 6 


Re / S v2 

r& = — ( ) 

4 x re f 


© Decoupling Requirement 


Ts » T L 


© Lower Bound of /? 




CRITERIA FOR 0 : LOWER BOUND 


0 Physical Requirement for Incompressibility 

The pressure wave must propagate much faster than the spreading of 
effect : 


Ti « Ts 

tl = time for pressure field to map the flow field 
ts — time for viscous effect to spread 


© Lower Bound of /3 


0 » 
0 » 



(Laminar) 

(Turbulent) 


CRITERIA FOR 0 : UPPER BOUND 


© Governing Equation in ‘Delta Form * 

(i + ~ r 0 + s v( A > n ~ r >) 

+^ c (A 3 b - r 3 )]) (Q n+l - Q n ) = RHS 

© Approximately Factored Form 

/ + £ - ro] • [/ + ^J n+l s n (A 2 n - r 2 ) 

. / + £ J»+^ c (A 3 n - r 3 )j A Q = RHS 


Higher order cross differencing terms are added. 



CRITERIA FOR 0 : UPPER BOUND 


O Added terms must be kept smaller than the original terms 

o Q JMi) > o 


O Upper Bound of 


(3h < 0(1) 




Re=1000, p=0.1 

(Recommended Range: 0.12< P < 10.0 ) 



Re=1 000, p=5.0 

(Recommended Range: 0.12< P < 10.0 ) 







LOGIRMSOO) 


Convergence History for Channel Flow 
Re=1000, Channel Length = 15.0 
(Recommended Range: 0.12< p < 10.0 ) 




Artificial Compressibility on Convergence 
Impulsively started circular cylinder to steady-state at Re=40 
(Recommended Range for p : 0.1< p<10.0) 

Stream Function Vorticity (Mehta) 






NUMERICAL ALGORITHM 


© LU-SGS Scheme 

(Implicit Lower-Upper Symmetric- Gauss-Seidel scheme) 

• Start from an unfactored implicit scheme 

{/ + £ [* { i + 6„B + f<c] ) (D“ +1 - D") 

= -A i - E„) + S V (F - F v ) + £ C (G - G„) 

• LU-SGS implicit factorization scheme 

LtL d ~ l L u = RHS 

where 

L, = 1+ jOr £ i + + <-,B + + S- ( C + 

Lj = i+ ~(A + - A + + b + +c + - d~) 

& 

L u = I + +A + + B + + C + ) 



Downstream Boundary Conditions 


0 Velocities 

• Extrapolate to update (u n )i^ a 

• Calculate mass flux out 


TO out 



• Then, obtain weighted velocities to conserve 


TO^tii 


Downstream Boundary Conditions (cont’d) 

0 Pressure 

• ^-momentum equation at L = L max — 1 and t = t* 

[d r uf + 9^! + + [9^C3] n = 0 

where e, = ± [wU { + (£,)<?] - 5(Vfc • Vfc)£* 

• Pressure corresponding to mass weighted velocities 


= P W — T" [(wC/3)** — (tlft/s)"] 


+ -(vc-vc) 


( & W \ h ( n 
\dc) “Vac ; 


(Variation in velocities normal to C-direction is assumed to be negligible) 


Downstream Boundary Conditions (cont’d) 
O Pressure (cont’d) 

• Progressive correction method 



Ip — f P n da 

Jtxit 

A momentum- weighted pressure is then formed as 


Skewness and Metrics 



— *1 = constant 


C- 

u -* 

ds l ■ ds L = (tj x 2 + I7, 2 )(y ) 2 

• bj -<*•»- ( 6 * + fe 2 )(^) 2 

« •/ 


dSi • C?32 

(rf«a 2 rf«a 2 ) 1/2 


= COS0 


* ^ 

Skewness and Viscous Terms 

O AF Algorithm for Incomp Navier-Stokes Eqs 

/ + ±J6 ( (Ai - e,t + e^A*)] [...] (...) ( Q n+1 - Q n ) = ^5 
Jiff 5 = - Ar J - £„!)...] 

0 Viscous Term 

E,x = j (V& • VZjI m 6 tj Q + ...] 

For orthogonal grid, (and for LHS) 

e,i = T(V^V^J m ^Q + ...) 

J 



Implicit Scheme 


© Linearize the right-hand side, using first order upwind difference 
and approximate Jacobians, resulting in a system of equations of 
the form 

B[T, 0, 0, U , 0, 0, X , Y, Z, 0, 0, V, 0, 0, W]AD = R 
0 Solve iteratively with line relaxation 

B[X, Y, ZJAD I+1 =R- TAD| - UAD|JL m 

- VAD‘j^ k - WAD[ JiM 


© Solve efficiently by 

1) First form and store T, U, V, W } X, Y, Z, and R for all grid points 

2) Perform LU decomposition of left-hand side, overwriting X, Y, 
and Z 

3) Iteratively form right-hand side and solve 


Boundary Conditions 


0 Inflow and Outflow boundaries based on Method of Characteristics 
0 Inflow Boundary 

• One upstream traveling characteristic 

• Three flow components specified, either 
All velocity components, or 

Total pressure and flow direction 

0 Outflow Boundary 

• Three downstream traveling characteristics 

• Static pressure 



Upwind Differencing I 


© Because artificial compressibility is used, convective terms form a 
hyperbolic system of equations 

• Upwind differencing used to follow propagation of artificial 
characteristic waves 

• Provides dissipation to suppress oscillations caused by the non- 
linear convective terms 

• Produces a more diagonally dominant system of equations than 
central differencing 

• More expensive to compute, but convergence improves consid- 
erably 

0 Flux-difference splitting with Roe averaging is used 

© Implemented as 3rd or 5th order accurate throughout all the interior 
grid points - no limiting 


ST 


Upwind Differencing II 


dE [j^t+l/2 “ &i- 1/2] 

” A£ 

E%+ 1/2 = 2 [^(A+i) + E(Di) - &+1/2] 

0 provides dissipaton, function of the flux difference A E± 

A E± 1/2 = A ± (D)(A + 1 - Di) 

© Roe properties are satisfied for D = 0.5(Dj+i + Di) 

A ± = XiAfXf 1 


= \ (^i + \/ A i + 



TIME-ACCURATE FORMULATION 

A 

© Time- accurate formulation not ds straight forward - no ^ term 
© First discretize the time term in momentum equations using second 
order three-point backward-difference formula 


{ 8U 8V 

\d( + dr, + 9( ) 

3u n+1 _ 4 ti n + u n_1 
2At 


n+1 

= 0 




© Iteratively solve these equations by introducing a pseudo-time level 
and artificial compressibility 


~ /gn+l.m+1 

A 

A t K 


_ £n+l,m^ — — y9V • ti n+1>m+1 

o»n+l,m 

— U n +^ m ) — 


- 4 u n + u n ~ l 
~2A t 


Resulting Systems of Equations 


© Steady-state formulation 



© Time-accurate formulation 


Itr 

J 


+ 



A D = 


J^Tl+lyTTl 


—(3 D n+1 > m - 4 D n + D n 
2AV 


') 


© Same residual R for both formulations, formed using 

• Upwind differencing for convective fluxes 

• Central differencing for viscous fluxes 


Summary of Time-Accurate Formulation of Artif icial 
Compressibility Approach 

• Time accuracy is achieved by subiteration 

- Discretize the time term in momentum equations using second-order 
three-point backward- difference formula 

3q"'-4q- + q ^_ _ 

2A t 

- Introduce a pseudo-time level and artificial compressibility, 

- Iterate the equations in pseudo-time for each time step until 
incompressibility condition is satisfied. 


1 ^/ « 

Tx « -i 


■ ) = -(rhs)~ tJ ~* - 


3 ^- 4 ?-+^ 


• Code performance 

- Computing time : 50-120 ms/grid point/ iteration 

- Memory usage: Line -relaxation 45 words/ grid point 

GMftES-ILU(O) 220 words/grid point 


i 




Pressure Projection Methods 


• Pressure is used as a mapping parameter to maintain incompressibility, 
that is, to maintain divergence free velocity field 

=> This is accomplished via derived equation, i.e. Poisson equation 
for pressure 

• The solution is usually done in multi-step (fractional step method) 

• Basic idea of this approach will be illustrated by presenting 

- A AAC method 

Issues n pressure field computation 

- A form of SIMPLE method 

- A generic pressure projection method via fractional step approach 

• A 3-D generalized method will be presented in detail 


Marker- and- Cell (MAC) Method 

Harlow and Welch (1965) 

O Poisson Equation for Pressure 
(Taking divergence of moment um eqn) 


where 


_ d dui _ 

P dxi dt dxi 9 


hi = 


dUjUj dTjj 
dXj dXj 


O Solution Procedure 

• Momentum equations are solved for velocity 

• Poisson equation is solved for pressure requiring divergence 
free velocity field at the next time level 





MAC Method 


© Grid 


»o+i 



f 




- 

1 

*- • — 
p 

♦ 





[ I 





staggered regular 


© Original MAC method : staggered Cartesian grid 

• Conserve mass, momentum and kinetic energy 

• Avoid odd-even point decoupling 

=> In generalized coord these advantages become unclear 
© Poisson solver requires large computing time 


MAC Method 

© Important to get numerical divergence free velocity field 
© Poisson Solver 

Method 1 : Use an exact form of the Laplacian 

V 2 p = g' 

Fourier transform of this 

-k P = g 

where 

p = Fourier transform of p 
fc 2 = fc x 2 + ky 2 + k z 2 

2 7T . , 

k x , .. = = wave number inx — direction 

g' = finite difference approximation to g 

n = —N/2, ... 0 , 1 , ~(N/2 — 1 ) 



MAC Method 


G Poisson Solver 

Method 2 : Use difference form of second derivative 

,6 2 6 2 6 2 x 

'6x 2 + 6y 2 + t>z 2 ’ P 9 

Fourier transform of this 




where 


(ki) 2 = Fourier transform of 


s 2 


SxiSxi 

( i.e. for a five point centred differencing ) 

(k) 2 = ^2 I 15 “ 16cos(A*i) + cos(2A*<)] 


MAC Method 

® Poisson Solver 

Method 3 : Use divergence-gradient operator 
Finite difference form of the governing equations are 


Sui 

It 



Sp_ 

Sxi 


= Gp 


Sui 

Sxi 


= Du{ = 0 


fajUj STjj 
fix jf Sxj 


where 



MAC Method 


0 Poisson Solver 

Method 3 : Use divergence-gradient operator (cont’d) 
Applying the divergence operator, D, 


DGp = 


6Dui 

6t 


+ Dh'i 



Fourier transform of this : —k'ik'ip = g' 
where (&j) 2 = Fourier transform of DG operator 
(i.e. for a five point central differencing) 


{hi) 2 = ^2 [65-16cos(AJfei) - 64cos(2AAr<) 
+16cos(2Afci) — coa(4AA:i)] 


MAC Method 


0 Poisson Solver 

• Fourier transform of momentum equation 



St 



To satisfy the continuity equation in grid space 



MAC Method 


O Poisson Solver 


• Substituting p fromm the three methods into the above 
For Method 1: 


l(k[u,) = (h[ - / 0 


For Method 2: 


L/t, f 

(K - * o 


For Method 3: 

Tt 


(^- k 3m= o 


Comparison of Laplacian Operator 






Semi-Implicit Method for Pressure-Linke d Equations 
Caretto, et al. (1972), Patankar and Spalding (1972) 


© Computing correct pressure for a divergence free velocity field 
in each step (i.e. Poisson solver for pressure) slows down the 
overall computational efficiency 


© For a steady-state solution, the correct pressure field is desired 
only when the solution is converged 

© Iteration procedure for the pressure am be simplified such that 
it requires only a few iteration at each tune step. 


=* The most well known of this approach is the SIMPLE method 


SIMPLE Method 


© Solution Procedure 

• Begin with a guessed pressure p* 

(Usually p n at the beginning or the cycle) 

• Solve momentum equation to get an intermediate velocity 

u** 

Ui* - Ui n = At[/cn(ti*,v n ,ti*,i>*) - ^-] 

i 

• Corrected pressure and velocity are obtained by 

P = P*+p' 

Ui — Ui*'+ Ui' 


Q 



SIMPLE Method 


© Solution Procedure (cont’d) 

• Relation between p' and v! 

- First, linearize the momentum equations 

- Then, drop all terms involving neighboring velocities 




(fen differencing scheme) 


S£ 

Sxi 


• Substituting these into the continuity equation, a pressure 
correction equation is obtained 


SIMPLE Method 


© Salient Features 

• This procedure in essence results in a simplified Poisson 
equation for pressure, which can be solved iteratively line- 
by-line. 

• The unique feature of this method comes from the simple 
way of estimating the velocity correction ttj. 

• This feature simplifies the computation but introduces em- 
piricism into the method. 

• Despite its empiricism, many computations have been done 
successfully using various forms of this method. 



Fractional Step Method 
Chorin (1968), Yanenko (1971), Marchuk (1975) 

© Time evolution is approximated by several steps. 

© Various operator splitting can be adopted by treating the mo- 
mentum equation as a combination of convection, pressure, and 
viscous terms. 

© The common application of this method is done by two steps. 

(1) Solve for an intermediate velocity field using a simplified 
momentum equation. 

(2) Compute pressure which maps auxiliary velocity onto a 
divergence-free velocity field. 


Fractional Step Method 


O Solution Procedure (an example) 

• Step 1: Calculate intermediate velocity, Hi 

(i.e. by a second order Adams-Bashforth method) 

~ *" - £ + + “*"> 


where 


Hi ~ 


Fractional Step Method 

O Solution Procedure (cont’d) 

• Step 2: Solve for the pressure correction 

^ n+1 - Uj _ 1 _ ^*) 

At 2 Sxi KV * J 


where 


_n in ^ 


This combined with continuity equation results in Poisson 
equation for pressure correction. 

VV + * = 


Fractional Step Method 


O Solution Procedure (cont’d) 


• New pressure and velocities are calculated as follows: 

p n+l =p n + W „+l _ r ) _ ^V 2 (^ +1 - r) 
Uj n+i = 


Fractional Step Method 


O Salient Features 

• Need special care for intermediate boundary conditions 
(Orszag et al., 1986) 

• As other pressure based methods, efficiency depends on the 
Poisson solver. 

• A multigrid acceleration is one possible avenue to enhance 
the computational efficiency. 



Pressure Projection in Generalized Coordinates 


• Approach in generalized coordinates 

- Finite volume discretization 

- Accurate treatment of geometric quantities 

- Dependent variables - pressure and volume fluxes 
Mass and momentum conservation 

- Implicit time integration 

- Fractional step procedure 

Solve auxiliary velocity field first, 

then enforce incompressibility condition by solving a Poisson 
equation for pressure. 

• INS3D- FS Code performance 

- Computing time : 80 ms/grid point/ iteration 

- Memory usage: 70 words/grid point 




Pressure Projection Method 


• Fractional-Step 

- Solve for the auxiliary velocity field, using implicit predictor step: 

+/»(«;> 

At 

- The velocity field at time level (n*l) is obtained by using a 
correction step: 

2-(u"-u,)=-Vp" +ft(0-V+M<0 
At 

- The incompressibility condition is enforced by using a Poisson equation 
for pressure [p* m P~* - P) 

V*p' = —V u 

At 




FORMULATION - II 


• Conservation of volume for time- varying cell 


8V 

dt 



dS • v = 0 


• Mass conservation 



dS-u = 0 


DISCRETIZATION - I 
Geometry 


• Closed cell 


£ * 

1 = faces 


= o 


• No gaps or 

overlapping cells 

Y, V = Vuui 

cells 




DISCRETIZATION - II 
Geometry (cont.) 

Conservation of volume for time- varying cell 
can be satisfied by 



DISCRETIZATION - III 
Mass Conservation 


(# • *)n - (3* • u) s + (S' 1 
• u)m — (S^ - u)p = 0 

Ut = S*-u 


u)e — (S 1 * * «)w 



IP = S' 1 - u 
U< = • tf 


DISCRETIZATION - III 
Momentum Conservation 




U* - momentum: 



since 


S* • u = U* 



DISCRETIZATION - IV 
Summary of properties 

• Dependent variables: U l , P 

• Staggered grid 

• Conservative scheme 

• Second-order in space 

• First-order in time 

• Implicit scheme 

• Fourth-order conservative numerical dissipation 


Test Problems for Code Validation 


• Two problems are selected to validate INS3D codes and compare 
artificial compressibility method vs pressure projection method 

- Wake and wing tip vortex propagation 

- Impulsively started flat plate 

Performance of any codes will depend on many factors beyond 
algorithm. However, it is assumed that both algorithms are reasonably 
well coded and optimized at a similar level 
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Vorticity- Velocity Method 
Fasel (1976), Dennis et al. (1979), ... , Hafez (1988) 

© Momentum equation is replaced by vorticity transport equation 

Bui du>iUj i dui 2 

dt dxj i 3 dxj 

where the vorticity u> is defined as 

u) — V x u 


© Taking curl of the above and using continuity equation, 

V 2 u = —V x a; 


© These are solved for velocity and vorticity. 

=>• Computational efficiency of this approach depends on the Poi- 
son solver. Overall performance in general three-dimensioned 
applications r emains to be seen. 



Engineering Applications 


• Applications 

- Viscous incompressible flows are encountered in many applications 

- Computational tools offer greater flexibility for resolving engineering 
problems, i.e. developing/designing flow devices 

- Issues in engineering applications are illustrated by presenting the 
following two real-world applications problems 

• Turbopump 

- Most of the material was present at 

the First MIT Conference on Computational Fluid and Solid Mechanics, Cambridge 
MA, June 12-14, 2001, entitled as •High-End Computing for Incompressible Flows* 

• Biofluid 

- Most of the material was present at 

the Sixth U.S. National Congress on Computational Mechanics, Dearborn, MI 
August 1-4, 2001, entitled as ‘Computational Hemodynamics Involving Mechanical 
Devices" 




Applications Point of View 


• Applications to Real-World Problems 

- N-S solution of full configuration was a big goal in the 80s 

- Numerical procedures and computing hardware have been 
advanced enabling simulation of complex configurations 

• Some Examples of Successful Applications 

- Components of liquid rocket engine 

- Hydrodynamics (Submarines, propellers, ...) 

- Ground vehicles (automobile aerodynamics, internal flows...) 

- Biofluid problems (artificial heart, lung, ...) 

- Some Earth Science problems 

• Current Challenges 

- For integrated systems analysis, computing requirement is very large 
=> Analysis part is still limited to low fidelity approach 

- For high-fidelity analysis, especially involving unsteady flow, long turn-around 
time is often a bottle neck => Acceleration of solution time is required 







Artificial Compressibility Method (INS3D-UP) 


• Time accuracy is achieved by subiteration 

- Discretize the time term in momentum equations using second-order 
three- point backward- difference formula 

v ,l -v±il = - (r ^ r 

2 A/ 

- Introduce a pseudo-time level and artificial compressibility, 

- Iterate the equations in pseudo- time for each time step until 
incompressibility condition is satisfied. 


p~' M ) = - fa" 

AT 


T -(q~'~'-q~ t ') = -(rhsr 
At 


3 -4 q* +q~ 
2 At 


> Code performance 

- Computing time : 50-120 ms/grid point/iteration 

- Memory usage: Line -relaxation 45 words/grid point 

GMRES-ILU(O) 220 words/grid point 


Pressure Projection Method(INS3b-FS) 


[VASA] 


• Approach in generalized coordinates 

- Finite volume discretization 

- Accurate treatment of geometric quantities 

- Dependent variables - pressure and volume fluxes 

- Implicit time integration 

- Fractional step procedure 

Solve auxiliary velocity field f irst, 

then enforce incompressibility condition by solving a Poisson 
equation for pressure. 

• Code performance 

- Computing time : 80 ms/grid point/ iteration 

- Memory usage: 70 words/grid point 





History of INS3D Development 


• Code 

1982-1987 Original version 
(Kwok, Chang) 


1988-1997 IN530-UP 

(Rogers, Kiris, Kwak) 
INS3D-LU 

(Yoon, Kwok) 

INS30-FS 

(Rosenfeld, Kiris, Kwak) 


• Applications 


SSME Phase U+ 
HGM Redesign 

* 

Inducer 


Flight engines 
1995 


Human implantation, 1998 
Oe Bakey/N ASA VAD 


1 99 8-Pres«nt Combined version 
(Kiris, Kwak) 


Advanced liquid 
sub -systems 


Current Challenges 


• Challenges where improvements are needed 

- Time-integration scheme, convergence 

- Moving grid system, zonal connectivity 

- Parallel coding and scalability 

• As the computing resources changed to parallel and distributed 
platforms, computer science aspects become important such as 

- Scalability (algorithmic A implementation) 

- Portability, transparent coding etc. 

• Computing resources 

- "Grid" computing will provide new computing resources for 
problem solving environment 

- High-fidelity flow analysis is likely to be performed using "super 
node" which is largely based on parallel architecture 




Supercomputing Memory and CPU Requirements 


,s r ' 

,0 F i 



THEORETICAL PROCESSOR SPEED, FLOPS 


IHEC TURB nPgL/TY 



Parallel Implementation of INS3D 



e INS3D-MPI 

(coarse grain) 

T. Faulkner A J. bodes 

• IN53D-MPI / Open MP 

MPI (coarse grain) ♦ OpenAV (fine grain) 
Implemented using CAPO/CAPT tools 
H. Jin A C. Kiris 


a 

M 


MPI 
• •a 

h 

Group 1 


Group 2 


Group N 



• IN53D-MLP 

C . Kiris 


Group 1 


Group N 










Tima (*#c) par harrton 


Parallel Implementation of INS3D 
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Initial start 


A 


Time Step 18: Impeller rotated 8-degrees at 1007* of design speed 




PRESSURE 


VELOCITY MAGNITUDE 



Summary of A Generic Turbopump Simulation 




Problem size: 

- 34.3 Million Points 

- 800 physical time steps in one rotation 
CPU requirement: 

- One physical time- step requires less than 20 
minutes wall time with 80 CPU's on Origin 2000. 

- One complete rotation requires one- week wall 
time with 80 CPUs. 

I/O: 

- Currently I/O is through one processor. Timing 
will be improved with parallel I/O since time- 
accurate computations are I/O intensive 
Parallel Efficiency: 

MLP/OpenMP version requires 19-257. less 
computer time than MPI/OpenMP version. 

Pin- to- node for MLP version reduces computer 
time by 407* 




Discussion on Numerical Procedures 


• Finite Difference 

- Based on Taylor series expansion => Requires smooth grid 

- Need special care for grid singularity 

- Generally easier to use fine grids near wall at high Reynolds number 

• Finite Volume 

- Formulation is more physical (conservation Of properties) 

- Viscous flux calculation is not as straightforward 

- Difficult to implement higher order schemes 

• In actual implementation, however, 

- These differences become unclear 

i.e. - FV in curvilinear coordinates requires tots of averaging depending on 
definition of variables such as staggered vs cell vertex arrangement 
- Both FD and FV implementations are very similar near grid singularities 

- Major dif f erences come from time integration scheme which also af fects 
the computational eff iciencies, especially, for unsteady flow computations 




Discussion on Applications 

am 

• Rapid turn around can be accomplished through the use of 

- Algorithm : convergence acceleration such as multi-grid, and GMRE5 

- Parallel implementation 

• Total process time can be reduced by 

- Automatic solution process including CAD to grid procedure 

• Need f urther development of methodology as well as physics modeling for 

- Deep understanding of flow physics such as unsteady flow characterization for 
better aeroacoustics modeling,and flow induced vibration 

=> Is LES method mature enough for this? 

- Need to matrix IT tools to flow simulation for smart flow control and optimization 

• Efficient extraction of information is still a challenge 

On top of all these we still need trained CFDers to solve many unsolved 
real world problems, for development of flow devices and for better 
understanding of flow physics 






Applications to Biofluid Devices 

9 

• Motivation 


• Mechanical Heart Assist Devices 

Computational Issues and Requirements 
Pulsitile Device 
Axial Flow Pump 


• Computational Approach for VAD Development 
CFD Technology Developed for Space Shuttle 
Design Improvements Using CFD 


• Summary and Discussion 




Mechanical Assist Devices 



• Motivation 

- Over 5 million Americans and 15 million people worldwide suffer 
from Congestive Heart Failure (CHF) 

- CHF patients are still treated with drug therapy, 
however, at late stage heart transplantation is traditionally 
the only treatment hope 

- Mechanical heart assist devices are being used as a temporary 
support to sick ventricle and valves as a 


" BRIDGE-T O-TR ANSPLANT" or "BRIDGE-TO-RECOVERY" 


- Need for assist devices is very high 
Permanent VAD need : 25,000-60,000 / YR 

Current valve replacement : 120,000/ YR 
Donor hearts available : 2,000-2,500 / YR 












DeBakey VAD™ Pump Schematic 


FLOW STRAIGHTENER 
(FRONT BEARING SUPPORT) 


fNDUCER/IMPELL ER 


MAGNETS EMBEDDED IN 
IMPELLER BLADES 


DIFFUSER 

(REAR BEARING SUPPORT) 


VAD BEARINGS 
ZIRCON1A SHAFTS IN OLIVE 
RINGS, AGAINST SAPPHIRE ENDSTONES 


MOTOR STATOR 


Issues in Axial flow VAD 


Problems Related to Fluid Dynamics 

- Small size requires high rotational speed 
Highly efficient pump design required 

-High shear regions in the pump may cause excessive blood 
cell damage 

Minimize high shear regions 

-Local regions of recirculation may cause blood clotting 
Good wall washing necessary 

>5mall size and delicate operating conditions make it difficult 
to quantify the flow characteristics experimentally 











DeBakey VAD Development Timeline 


• CFD Assisted Design 

1994 -Kiris and Kwak continued design changes 

=* IMPROVE BEARING, HUB AND HUB EXTENSION DESIGN TO REDUCE BLOOD CLOTTING 






DeBakey VAD Development Timeline 



• Animal Tests 

1995 - Animal implantation: passed two-week requirements 

1996 - Full! design rights are granted to MicroMed, Inc. to produce the pump 

Began using bio- compatible titanium replacing polycarbonate 

1997 - Configuration design finalized 
















DeBakey \MD 







VIE 201 & 202 at lunch 



• Human Implantation in USA 


1999 - US Patent is granted for the device on September 9, 1999 

2000 - Over 30 patients have received the device 

The longest successful trial period to date in human is approaching 
1 year (as of July 2001) 


2 







NASA/DeBakey VAD 
Accomplishments to date (7/1/01) 



A patient In Munich fully mobile and 
discharged awaiting transplant 



A patient I n Zurich, owner/chef of a 
restaurant worked dally until hfs 
transplant 



90+ patients implanted 

Number of patients currently ongoing with device 
(longest patient supported with the device is now 
approaching one year) 

U5 trial 

Approved for 20 patients in a multi-center trial 
European trial 

Received *CE mark' (the EU equivalent to FDA approval) 
Results to date 

Favorable compared to existing VADs 

Small incidence of thrombus is being investigated 

=> Further computational support is essential 



The first patient 
In Houston with 
Ors. OeBakey 
and Noon on her 
discharge day 
after t ra replant 



Summary and Discussion 



Computational approach provides 

- a possibility of quantifying the flow characteristics: especially valuable for 
analyzing compact design with highly sensitive operating conditions 

- a tool for conceptual design and for design optimization 

* CFD + rocket engine technology has been applied 

- to modify the design of NASA/DeBakey VAD which enabled human implantaion 

* Computing requirement is still large 

- Unsteady analysis of the entire system from natural heart to aorta involves 
625 revolutions of the impeller 

- During on heart beat, impeller has 125 revolutions 

- With 1024 processors of Origin, one simulation (with several heart beat) from 
heart to aorta can be completed in one month 

* Further study is needed 

- to assess long term impact of mechanical VAD on human body, which requires 
modeling flexible wall and non-Newtonian effect among other things 

■ There exist some gaps between 

- CFD (assuming TT is a part of CFD applications) and biomedical expertise 





^ Progress to Date 

w tan* tar 

• CFD has pioneered the field of flow simulation for 

- Obtaining engineering solutions involving complex configurations 

- Understanding physics (critical to mission success) 

• CF D has progressed as computing power has increased 

- Numerical methods have been advanced as CPU and memory increase 

- N-S solution of full configuration was a big goal in the 80s 

- Complex conf igurations are routinely computed now 

- DNS/LES are used to study turbulence, 

(but not to resolve mysteries of turbulence) 

• As the computing resources changed to parallel and distributed 
platf orms, computer science aspects become important such as 

- Scalability (algorithmic & implementation) 

- Portability, transparent coding etc 

- Coding paradigm is being changed (i.e. object-oriented program...) 


Simulation and Control Roadmap 


I \ AS\1 


CHIMERA TOOLS 
Overset connectivity 
•surface & volume grids 
•grid manipulation 


Post Processing 
-feature extractions 
•data compression 


m 





Examples of Current Capability 


• Algorithmic advances include 

- Discrete models : 

Various artif icial dissipation models 
Unified formulations, e.g. preconditioning 
Unstructured methodology 
Various gridding strategies 

- Solution methods: 

Explicit/Implicit 
Preconditioning, dual-time 
Multi-grid 

• Successful application of CFD to engineering problems 

- High- lift configurations 

- Multiple bodies in relative motion 

- Components of propulsion system (both aero 4 space) 

- Maneuvering vehicle 

- List goes on 


• EXAMPLE LANDING CONFIGURATION 

- 2 2. 4M mesh points 

- 79 zones 

- 201 C90 hours for convergence 
(Lift within 2% of experiment) 

=> Small geometric variations have a major impact, 
particularly near maximum lift 
=> Grid density study was performed 
=> Accuracy of physical modeling needs further assessment 




Stuart Ro?**, NASA Amaa - AST/IWO High Lift 



Overset Technology for Complex Configuration 




• Overset (Chimera) Grid Approach 

NASA Ames Developed CFD Tools 

- OVERFLOW Navier-Stokes Flow Solver 

- Chimera Grid Tools: Pre- and Post 
Processing 

- Enabling flow simulation technology for 
complex configuration and unsteady 
flow involving bodies in relative motion 


OVERFLOW+CGT: 1996 NASA Sdtware Of th# 

Year Hononbto Mention 

PL0T30 Visualization Software 

1993 NASA Spaca Act Award 

FAST Visualization Softwara 

1995 NASA Software of the Year Award 


Overset Technology for Complex Conf iguration 



X38 B52 
Drop 


• Wide Range of Applications 

- Spacecraft ascent and descent (still need 
better chemistry) 

- Propulsion (limited scope) 

- Aircraft 

- Hydrodynamics (limited unsteady capability) 

• Current Development 

• Working toward fully automated grid 
generation 

- Steady and Unsteady capabilities (in aM 
speeds) 

- Bodies in Motion, 6 DOF 

- Non-equilibrium chemistry 


-\ .. 


Turbopump 




777 




Scalability of Parallel Cart3D 



• Parallel Version of Cart3D 

• Parallelization through domain decomposition 
and explicit message passing. 

- Domain decomposition technique based upon 
space-filling curves permits domain 

’ decomposition to be performed in parallel 
and on-the-f ly (at runtime) 

- Parallel speed-ups of ~53 on 64 processors. 

- Combined with robust multigrid to offer 
exceptionally fast convergence across the 
range of Mach numbers. 

• Validation using Citation Twin-engine 

Business Jet 

- 1.42M cells. Mach 0.84 alpha 1.81deg 


Unsteady Simulation of A Generic Turbopump 


Overset grid set up 


• Status 

- 34.3 Million Points 

- 800 physical time steps / rotation 


1st time step 



Velocity vectors: 

Impeller started at 1 0 % actud 
rotational speed. 


18 th time step 


Velocity vectors: 

Impel tor rotated 8-degrees 

at 100% of design speed 


=>One physical time-step requires 
less then 20 minutes wall time with 
80 CPU’s on Origin 2000. One 
complete rotation requires one-week 
wall time with 80 CPUs. 

=> Currently I/O is through one 
processor. Timing will be improved 
with parallel I/O since time-accurate 
computations are I/O intensive. 

• Current Goal 

- Further automation of moving grid 

- Parallel I/O 

=> One rotation per day 



Parallel Implementation of Flow Solver 




INS3D-MLP / 40 Groups 


A Generic Turbo pump 
1 14 Zones / 34.3 M grid points 


g BOO 
S 400 
1*0 
* 


I 


-■—02000 no -pin 
— D— 02000 pi n 
— 03000 no -pin 
-□-03000 pin 


Per processor Mflop Is 
between 60-70. 

Code optimization for 
cache based platforms 
Is currently underway. 
Target Mflops is to 
reach 120 per 
processor. 

Increasing number of 
OpenMP threads Is also 
the main objective for 
this effort 


200 300 


Number ol CPUs 


Static/Dynamic Stress Analysis 
for Flow-Induced Vibration 


NEW CONDITIONS 



INTERFACE 


STRUCTURAL 


ZIPPER/GRID 

A. 


LOADS 


CFD GRID 
PRESSURE 
TEMPERATURE 


FLUIDS 

INS3D 


FEM GRID 

I 


STRUCTURES 

NASTRAN/ANSYS 


STRESS /FATIGUEl 
ANALYSIS I 


USE PSE/IPG TOOLS (GROWLER/SCIRUN/PUNCH) 



Performance (QFLOP/s) 


Examples of Current Capability: OVERFLOW-MIP Performance 


System: 512 CPU Lomax (300 MHz Origin 2000) 


Problem: 35 M Points, 16p Zones 


$teger (32^it) 


Lomax (j&4bit) 




« ✓ i 

if*** - > 

- " : 


Steger (64: 


- FV^e-CeveH-MHe stone - 


: 16 CPU C90 OVERFLOW limit 
, Origin limit fpr Standard fMaaaa 

256 

Number of CPUs 


Examples of Current Capability: OVERFLOW-MLP Performance 


• Origin 2000 (64 bit) performance is dramatically better than full C90 

• OVERFLOW 16CPUC90 = 4.6 SFLOP/s 

• OVERFLOW 256 CPU 02 K (250MHz) = 20.1 SFLOP/s 

• OVERFLOW 512 CPU 02K (250MHz) = 37.0 SFLOP/s (cluster) 

• OVERFLOW 512 CPU 02 K (300MHz) = 60.0 SFLOP/s 

• Performance/Cost Advantage of Steger/Lomax over C90 

• OVERFLOW = 256 CPUs are 4.4x faster @ 4.5x Cheaper = 23x 

• OVERFLOW = 512 CPUs are 13.0x faster @ 2.6x Cheaper = 33x 


• Performance gains for small changes in code 

• ~1000 lines of changes (<1% of total code) 



Status of CFD Research 

mrnmmtcmm 

“Can do it all* message was propagated in the past, but 
CFD did not replace Wind Tunnel => CFb was oversold! 

Of course, further research will create advances with across the board benefits; 

• Algorithm 

Convergence acceleration, Robustness, Error estimation 
Grid related issues, adaptive grids 

• Physical modeling issues 

Turbulence, Combustion, Multiphase, Cavitation, Spray, Plasma etc. 

• Solution Procedures 

Automation: C A b- Grid -Solution -Feature extraction 

• Applications 

Rapid turn around for complex conf igurations 

besign and product development - we still need trained CFDers 

=> Outsourcing makes sense 

However, sponsors are likely to view these as “incremental advances." 


Current Needs in Flow simulation 


Work environment is different now 


- Tremendous information is availabie 


- Single-handed code development is rapidly becoming outdated (CFD 

discipline as defined in the past is disappearing) 

- Problem solving environment is more collaborative 

=* Requires software engineering to mitigate risks: 

Legacy software handling tools 
Visualization 
Data base handling tools 
Object-oriented coding 






Technology Need 


Integrated solution for assessing the total system performance, life cycle 
and safety can very well be the next challenge 

e.g. Need a more complete picture of entire design space not just one design 
Some challenges specific to CFD are: 

- Physics-based simulation for more predictive capability 

- Integrated analysis 

e.g. multi-discipline, performance for entire flight envelope 

- IT tools can be used to integrate CFD, experiments and flight tests 
e.g. virtual flight 

=> Requires : Many simulations which will be put into data base, and 
data base management tools, query tools to extract desired info 

1 Validation is an issue 


Example: Real Gas Effect Model 


ELECTRONIC STATE FOR N2 



Excited State 


Ground State 


U U U M M M 


• Current Model 

Euler : 

- Do not require knowledge of internal 
molecular structures and intermolecular 
potentials 

Navier-Stokes : 

- Molecules are structureless 

- Transport properties are based on a 
single intermolecular potential 

- Collisions are assumed to be elastic 

Non- equilibrium flow equations: 

- EOS for each species is based on equil 
distributions over many internal states 

- Reaction rates account for ground states 

- Empirical intermolecular potential is used 




x 10-* 


Example: Real Gas Effect Model 


ELECTRONIC STATE FOR N ATOM 
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Excited State 


Ground State 


• Proposed Approach 

- Based on more accurate solution of 
known microscopic equations, develop 
better macroscopic equations 

* Derive micro eqs and constitutive eqs 
from Bloltzman eq (inelastic collision) 

- Obtain state- to-state rates and product- 
state distribution functions 

=>Provide macro properties to be used in 
CFD codes 

• Impact 

- The results are more accurate physics- 
based representation of macroscopic 
properties (from current curve fitting) 

- Applicable to high-speed planetary 
reentry 


Example of Next Level CFD Problem 


Typical Descent Trajectory 


II ! Ill ; IV ; 


Non- Continuum 
Region 

Continuum 
Region i 

* ■ I 


V VI • Tough Problems: 

£ Physics-Based Scientif ic 

TrajepVy Computing + CFD 


Physical Modeling 


Aerodynamics* I Aerothermodynemics 

- PwtidGM** * Equab RtadtogGaa ► 

^ * J **-Equ*) >*- 

m • m \ j ! \ * Radiation 

4 0 12 16 20 24 

M - Mach Number 

Generic Vehicle 


• Big Impact on Engineering : 

Major contribution in developing a 
vehicle/system 




Examples of Target Problems 


• There are many unsolved challenges in developing flow devices 

- Compressor rotational stall 

- Turbopump system in rocket engine 

- Jet noise 

- Maximum lift of high-lift system 

- Rotor-based propulsion system 


Examples of Target Problems 


• There are a wide range of challenging applications in non-aerospace 
Earth simulation 

- Climate prediction 

- Local / regional model 
Biofluid 

- Flow-related problems in human body (e.g. heart, lung, hemodynamics): 
fundamental understanding / biomedical practics 

Engineering and product development 

- Automobile 

- Naval hydrodynamics 

- Chemical engineering 

- Manufacturing processes 



Strategy 




• Bottom line is "money" 

- Traditional "CFD* research we used to know is probably over. 

• Short term 

- Very often, need to satisfy "paying customer" 

- Software for any conventional engineering is most likely available 

- Risk might come from lack of engineering knowledge by CFDers 

• Long term 

- Basic research is an investment for the future 

- However, need a "vision* for planning 
Need to identify high payoff areas 




